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Abstract 

It is shown that the vibrational wave packet relaxation of initially co- 
herent (displaced) states as well as the quantum superposition of coherent 
states in heat bathes with the different spectral densities exhibit a number 
of peculiarities compared with the cases of linear phase-sensitive relaxation, 
quadratic relaxation, and relaxational dynamics described with the use of 
Zurek's environmentally-induced pointer basis. A strong dependence of the 
relaxation rate on the position of the spectral density maximum of the bath is 
found. The differences can be used for the discrimination of the mechanisms 
of the molecule-environment interactions. 
PACS number(s): 82.29.Rp, 33.80.Be, 32.80.Bx, 31.70.Hq 
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I. INTRODUCTION 



The behavior of many quantum systems strongly depends on their interaction with the 
environment. It is important to take this interaction into account to reahstically describe 
systems hke, e.g., vibrational levels in a big molecule, the quantized mode of an electromag- 
netic field, or a trapped ion. The rapid development of experimental techniques in these 
and other branches of physics and chemistry leads to an increased interest in theoretical de- 
scriptions of possible experiments by numerical calculations. Additionally our calculations, 
which concern such systems, allow us to regard the still existing question about the border 
between classical and quantum effects from a new point of view. For the systems, which are 
described here, the border to disappearance of quantum effects has been estimated. 

One of the fundamental questions of quantum physics is to understand why the general 
principle of superposition works very well in microscopic physics but leads to paradox situa- 
tions in macroscopic physics such as the Schrodinger cat paradox [|l| where the cat can exist 
in a superposition of the states dead and alive. One possible explanation of the paradox and 
the nonobserving of a macroscopic superposition is that |^ systems are never completely 
isolated but interact with an environment, that contains a large number of degrees of free- 
dom. Interactions with the environment lead to continuous loss of coherence and drive the 
system from a superposition into a statistical classical mixture. 

The interest in the decoherence problem is explained not only by its relation to the 
fundamental question: "Where is the borderline between the macroscopic world of classical 
physics and microscopic phenomena ruled by quantum mechanics?" , but also by the increas- 
ing significance of potential practical applications of quantum mechanics, such as quantum 
computation and cryptography [3, 4]. 

There are a number of propositions how to create the superposition states in mesoscopic 
systems, or systems that have both macroscopic and microscopic features. Representative 
examples are the superposition of two coherent states of an harmonic oscillator 

= A^-i + e^<^|-a)) (1) 
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for a relatively large amplitude {a ~ 3 -^ 5). Here, \a) is a coherent state and = 



recently for the intracavity microwave field |^ and for motional states of a trapped ion 0. 
Additionally, it has been predicted that superpositions of coherent states of molecular vi- 
brations could be prepared by appropriately exciting a molecule with two short laser pulses 
and the practical possibilities of realizing such an experiment have been discussed P|. 
In this scheme the quantum interference would survive on a picosecond time scale, which is 
characteristic for molecular vibrations. 

From the theoretical point of view, quantum decoherence has been studied extensively 
[2, 9-16]. Most efforts focused on the decoherence of the harmonic oscillator states due 
to the coupling to a heat bath, consisting of a large number of oscillators representing 
the environment. The system is usually described on the basis of the master equation 
for a reduced density operator. There are two general approaches for this method. One 
adopts the Markov approximation together with the rotating wave approximation (RWA). 
That means that all details of the complex system-environment interactions are neglected 
and relaxation is described by characteristic decay constants. The system-bath interaction 
is phase- insensitive because of the RWA. In another approach, according to Zurek 0, the 
coupling with the environment singles out a preferred set of states, called "the pointer basis" . 
Only vectors of this basis survive the quantum dynamics. The vectors of the pointer basis 
are the eigenvectors of operators, which commute with the (full) interaction Hamiltonian 
of the system. This basis depends on the form of the coupling. Very often this pointer 
basis consists of the eigenstates of the coordinate operator. The density operator describing 
the system evolves to diagonal form in the pointer basis, which is usually connected to the 
disappearance of quantum interference. The two approaches give different pictures of the 
same decoherence processes. 

One of the goals of this contribution is to present a consistent analysis of the deco- 
herence on the basis of a density matrix approach starting from von Neumann's equation 
for the density matrix of the whole system, i.e. the microscopic quantum system and the 




These states have been observed 
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" macroscopic" environment . 



II. GENERALIZED MASTER EQUATION 

Let us consider a single molecule vibrating in a one-dimensional harmonic potential. The 
molecule interacts with a number of harmonic oscillators modeling the environment. In the 
interaction Hamiltonian 

//es = ^ E i^t + h) (« + (2) 

a (a"*") are annihilation (creation) operators of molecular vibrations with the frequency cu, 
(b^) operators for the environmental vibrations having the frequencies uj^. is the coupling 
between them. Starting from the von Neumann equation 

P--\[H.P\. (3) 

where the Hamiltonian 

H^Hs + He + Hes (4) 

contains the molecular system Hs — hcu {a'^a + 1/2) and the environment He — 
Eg hiu^ (6+6^ + 1/2), one can formally rewrite the equation for the reduced density matrix 
a = trE (p), which is averaged over the environmental states [17, 18] 

a = -1 [Hs, a] + e-^^«*/'^I)(t, 0) {D{t, 0))"' e^^«*/V, (5) 

where 

L'(i,0)=tr(f'exp^*cir(^-^L(r)) •pb(O)) (6) 

is an evolution operator, averaged over the initial states of the environment. 

L(r) = [i^Es(r), ] (7) 

is the Liouville operator in the interaction representation and T is the operator of chrono- 
logical ordering. Supposing that the initial states of the bath oscillators are thermalized 



(p^(0) ~ exp{—hu!^b^b^/kBT)), and restricting a cumulant expansion of D to the second 
order cumulant [17, 18], we obtain the non-Markovian master equation 



a — —lu 



a'^a, a 



+ Rcr, 



(8) 



where the action of the relaxation operator R is defined by 



Ra = 



a"^ + aj + j^a"*" + a,a (^An + 74^_|_;^^ 



(9) 



Here, the operators An and An+i are defined by linear combinations of the operators a and 

a"*" as 



An ^1n{t)a + %{t) a+. 



with the functions 



-fn+l{t) = f Rn+l{T)dT = ^ + 1)' 



(10) 
(11) 



(12) 



7„(t) = / Rn{r)dT = Y.Kln^- 



-i{tjj^—u))t -j^ 



(13) 



7n+i(t) = / i?„+i(r)e-2--dr = + 1) — 



g-i(a;{+a;)t _ ^ 



(14) 



7> 



.(t)= / i?„(r)e-2--ciT = ^/r|n^— . 
io ^ -z(a;^ + a;) 



(15) 



Note, that and i?n+i are the correlation functions of the environmental perturbations 
-Rn(i) — -f^l'T'^ exp (— i (o;^ — lo)t\ where denotes the number of quanta in the bath 
mode with the frequency a;^. 

The functions 7 correspond to the friction coefficient in the classical limit. The first 
two coefficients 7^ and 7^+1 strongly depend on the coupling constant for frequencies 
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uo^ c:^ uo and on the number of quanta in the bath modes with the same frequencies, whilst 
the coefficients 7„ and 7„+i are very small for all frequencies. 

The obtained master equation describes different stages of vibrational relaxation. 
The initial stage is defined by a period of time smaller than the correlation time Tc of the 
environmental perturbations. This time can roughly be estimated from the width Auj of the 
perturbation spectrum K"^, for Tc ~ 1/ Auo. For such small times one can write the master 
equation in the following form 



-lOJ 



a a, a 



+ r { (a+ + a) cr, a+ + 



+ 



a + a^a [a 



(16) 



where F = (2ri^ + 1) is a real constant. As follows from Eq. (|T6|), the pointer basis 

for this step of relaxation is defined by the eigenstates of the position operator Q ~ a'*' + a. 
Another period of time, for which the form of the relaxation operator R according to Eq. 
is universal, is the kinetic stage, where t ^ Tc and the Markov approximation becomes 
applicable. In this stage the master equation has the form 



cr 



-tLU 



a a, a 



(17) 



+7 I (j^n + l)a + na^^ a,a'^ + a + + a,a (j^n + l)a + na"*"^ | 



where 7 = irK'^g is the decay rate of the vibrational amplitude. Here n = n^, K = 
and the density of bath states g = g{uj^) are evaluated at the frequency u = of the 
selected oscillator. It should be stressed that Eq. ([T7| ) differs from the usual master equation 
for a damped harmonic oscillator for derivation of which the RWA is applied [5-8, 10-13]. 
This phase- sensitive relaxation leads to new effects: classical squeezing and a decrease of 
the effective harmonic oscillator frequency |1^. In between there is a time interval, where 
relaxation is specific and depends on the particular spectrum of K'^. 



III. ANALYTICAL SOLUTION FOR WAVE PACKET DYNAMICS 



The solution of the equation of motion of the reduced density matrix can be conveniently 
found using the characteristic function formalism, introduced in [|l9l, which enables us to 
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use the differential operators ^ and ^ instead of a"*" and a. Multiplying both sides of 
Eq.(|D with the factor / = exp (Aa"*") exp (— A*a) and taking the trace one can rearrange all 
terms into such a form that a and a'^ precede the appropriate exponent. For this opera- 
tion we change the order of operators using the expression a exp (Aa+) = exp (Aa+) (a + A) 
to make the transformation a+exp(Aa''") = ^exp(Aa+). After that every term can 
be represented by the normally ordered characteristic function F = tr(/cr) upon which 
one of the differential operators acts. We obtain, e.g., ti {[a~^a, a] f) = (A*^ — A^^ F, 
tr ([a+(j, a] f) = A* (^-^ — A*^ F, and tr ([aa, a^] f) = — A*^F. Such manipulations lead us 
finally to the com^/ex-valued partial differential equation: 

F = - [^u;X* + /. (A + A*)] ^F + [zuX - /x* (A + A*)] ^F 

- (A + A*) (z/A* + z/*A) F, (18) 

where 

Ki)=7:W+7n+i(t) (19) 

are relaxation functions. We can solve Eq. (^) by using the integral representation for the 
characteristic function F (A, A*, 0) = tr (^e"'^^e'°'^*^ which formally allows us to describe the 
nondiagonal density matrix. Below the notation / / dadPc{a, P) = {a \P) is adopted. An 
initial characteristic function 

F(A,A*,0)=y J dad(3e"^-f^^' c{a, (3) (20) 

will evolve in accordance with Eq.(p!8|) as 

F(A,A*,t) = / / dadpexp irt;^)(t)A'" (-A*)") c{a,f3). (21) 

\m,n ) 

We restrict the cumulant expansion to the second order, i.e. m + n < 2 . For a wide class 
of initial states (coherent, thermal, squeezed, etc.) high order cumulants vanish and our 
approximation becomes exact. The cumulants could hold nondiagonal information, such as 
the density matrix, in relevant cases we stress it with the upper index [a] or (/5). 



The functions Kmnif) in Eq.(pTD are given by the solutions of the sets of equations 



K 



(/3) 
01 

(a) 
10 



(22) 



= 2Rez/ - 2 (Re/x) Ku + 2fxKo2 + 2fi*K2o, 
K20 = -y* + ^lKll + 2{iuj- 12*) K20, 
K02 = -V + /i*i^ii -2{iu- /i) Kq2, 



(23) 



(24) 



with the initial values 

ir£)(t = 0)=a, 

ir(f(t = 0)=/5, 
Ku (t = 0) = 7^20 (t = 0) = 7^02 (t = 0) = 0. 
For the special case (3 = a*, these initial conditions represent a coherent state with amplitude 
a. This solution can be used for the construction of wave packets in different representations. 
Here, we will discuss the coordinate representation, in particular the dependence of the 
probability density P on the vibrational coordinate Q and on time 



where 



1 r'^ 

ZTT J-oo 



X (A, t) = tr (^e'^('^^+'^)a(t)) = e'^'/'F {tX, -tX, t) 



(25) 



(26) 



is a characteristic function for the position operator. Evaluating the integral ( p5D we finally 
obtain 



P{Q,t) = J J dadpP'-'''^^ {Q,t) c{a,(3), 



(27) 



where 



P("'^) {Q,t) 



2^W^) 



exp 



Q - g("'^) (t) 



(28) 



and 

Q(°'/^) it) = K^^ it) + Ki^ it) , ^^^^ 
V{t) = l + Kn {t) + K2o {t) + Ko2 (t) . 

For the case P = a* the function Q*-"'"*-* denotes the expectation values of the coordinate 
operator of the coherent state, V is the broadening of the Gaussian packet (^). The 
distribution P can be used for the investigation of relaxation dynamics for any initial state 
of molecular vibration, but it is best suited for studying the evolution of states prepared as a 
superposition of coherent states. Below, we will discuss the relaxation dynamics of initially 
coherent states and the superposition (|l|) of two coherent states. 

IV. RELAXATION DYNAMICS OF COHERENT STATES AND THEIR 

SUPERPOSITIONS 

A. Coherent states 

Because of the importance of this case we will use it to investigate the relaxation dy- 
namics described by the master equation @. For coherently exited states |ao) the initial 
characteristic function is F (A, A*, 0) = exp {a^X — ao^*). Thus, the initial values for Kiq^ (t) 
and i^if (t) are and Oq- In the first stage of relaxation, when t <^ t^, the relaxation func- 
tions are fi (t) = 0, u (t) = Tt. Therefore, the solution of the system of Eqs. p^ , p3[) 
gives 

V ; ^^^^ 

Even in this early stage there is a small quadratic broadening of the wave packet 
p("o."o) (^Q^t'^j without changing its mean amplitude. After the intermediate stage of re- 
laxation the solution of the system goes into the Markovian stage of relaxation, where the 
master equation ( p!7D works. For this stage the solution of Eqs. (p2| , p3D reads 
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2Re {aoz [t)) e-^\ 
- + n-ne~^^^ 1 + 



V{t) 





(1 — cos 2Ljt) + — sin 2ut , 



(31) 



where 



= cos cut + {■y/oj) sin cut + i{u/u) sinut, 



UJ = Juj^-j2^ 



(32) 



n = 



[exp{huj/kT) - 1] 



Equation (|32|) demonstrates the decrease of the effective harmonic oscillator frequency due 
to the phase-dependent interaction with the bath. The broadening of the wave packet V{t) 
can be also considered as uncertainty of the coordinate. We would like to discuss both its 
increase and oscillation, described by Eq.(31) and displayed in Fig. 1. The increase of V(t) 
appears due to absorbtion of quanta from the heat bath. Such an effect is rather clear 
and can be obtained already using RWA. The oscillation of the broadening deserves more 
intent attention. Some authors [21, 22] have predicted the oscillation of second moments, 
using Green functions formalism. One feature of prediction [22] is that it assumes initially 
a squeezed state. Our prediction, represented in explicit form in Eq. (31) goes further. We 
would like to underline that in our case such an effect appears even for the usual coherent 
state as the initial one. In spite of the fact that this oscillation does not present quantum 
squeezing, because the width is never smaller than the ground state width, we stress the 
derivation of wave packet broadening oscillation in Eq. (31) as one important achievement 
of the present paper. 



For the initial superposition of two coherent states (|T]) the normally ordered characteristic 
function consists of four terms: 



B. Superposition of two coherent states 



F(A,A*,0) = 



a* ,a -1 —a*,— a \^ ^ a*, 



—a 



) 



(33) 
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with Fa^p = e^^ ^'^* . The first terms describe the mixture of two coherent states, |a) (a| + 
|— a) (— «|. The last two terms correspond to the quantum interference, i.e., they refiect the 



coherent properties of the superposition. In accordance with Eqs.(p^)-p9|) and the equation 
of motion (0) of the reduced density matrix this initial state evolves as 



P (Q, t) = — [p("*'") (g, t) + p(-"*'-") (g, t)] + p,„, (g, t) 



with the interference term 



Pint (g, t) 



-2|q|' 



e*^p("*'-") (g,t) + e-*^p(-°*'") (g,t) 



where P^"'^) (g,t) is given by Eq. (p8D with 



g("*.-) (t) = -g 



(-a*, -a) 



(t) = 2'Re{az{t))e 



-7i 



(34) 



(35) 



(36) 



g 



(t) = -g(-"*'°) (t) = 2ilm (az (t)) e 



-7* 



(37) 



(t) and z{t) are defined by Eqs.(pT|, |32| ). Rewriting the real part of Eq. (|35|) we finally 
obtain 



Pint (g, t) 



1 1 / 2^4[Im(«^(t))]^e-27*-g2 
exp — 2 |q;| H ^ 



iV2 



4V^(t) 



(38) 



Figure 2 illustrates, how the superpositional state (1) evolves in time in accordance with 
Eq. (34). It follows from Eq. (38), that the interference term describing quantum coherence 
in the system is only significant when 



[Im(az {t))f 
V{t) 



2 lal 



(39) 



This is true for 7^ <^ 1, and moreover, when the two wave packets of the state (|I]) come 
close together (i.e., at the moment when z (tj) i). Expanding Eq. (|39|) into a series at these 
points and taking into account 7^^ -C 1, it is easy to see that 
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,2 , [lm{az{ti))f 



Pi„t ~exp(^-2|ar + ^ (1 _ 27t,)J . (40) 

The decoherence is due to two reasons: The spreading V (tj) = 1/2 + AV (tj) n of the wave 
packet due to thermal excitations by the bath, and amphtude decoherence. The second 
means, that even for the case n = and Im {az (ti)) = \a\ quantum interference disappears 
exponentially with the rate 

iJec^2|a|S (41) 

This result obtained by Zurek p[ is the main reason why quantum interference is difficult 
to observe in the mesoscopic and macroscopic world. For example, a physical system with 
mass Ig in a superposition state with a separation of 1cm shows a ratio of relaxation and 
decoherence time scales of 10'^''. Even if our measuring device is able to reflect the quantum 



properties of the microsystem, nevertheless objectiflcation ||2^ occurs due to the coupling 
between the meter and the environment. The fundamental result of Eq.(|41|) is obtained 
no matter which approach is used, e.g. the RWA or Zurek's pointer basis approach or 
the self-consistent description of the present paper. However, the time dependence of the 
superposition terms of the distribution of Eq.(38) differs a little bit, which can be seen in 
Fig. 3. As quantum interference is more sensitive, we have used it for comparison of three 
different approaches to the present problem. Figure 3 shows the difference between the time 
evolutions, which result from Eq. (0), from Eq . (|3lD , and from the corresponding result of 
the RWA approach. 

These differences arise due to the fact that during the relaxation there is no constant 
pointer basis for all steps of the evolution. As follows from Eqs. ([16], |1^), this basis changes 



from the position eigenstate basis in Eq. (ITB) to the more complicated basis in Eq. (|T^. 



C. Partial conservation of superposition 



As mentioned above the interference term (35) of the wave packet decays inevitably. 
However, there exist proposals, how to slow down this process, reducing the consequence of 
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the interaction. In particular, one can consider a situation, when the majority of the bath 
modes have an off-resonance frequency in respect to the selected system. In such a system 
the processes which are induced by the interaction with bath, namely loss of amplitude and 
phase must be delayed. The hope to observe the conservation of a superpositional state as 
a consequence of the mentioned properties is the reason to put attention on such a system. 
As a simplest example we take the quantum system surrounded by harmonic oscillators 
with twice the system frequency cu^ — 2a;. It provides processes, when the system loses 
2 quanta and the bath obtains one quantum. Although this is still a resonace situation, 
we demonstrate below that it leads to unusual behavior compared to the usual situation 
cu^ — cu, discussed above. Describing such a situation in the RWA we rewrite the interaction 
Hamiltonian (2) in the following form 

Hi^h^K^ (b+a' + k{a+f) . (42) 

Applying then the above mentioned formalism of the evolution operator (6) we obtain again 
the equation (8) for the reduced density matrix a of the selected system, but with some 
changes of the relaxational part (9) 

Ra = T{n^ + 1) { [aV, {a+f] + [{a+f, cra^] } 

+ Tnj { [(a+) V, a^] + [a^ a{a+f] } , (43) 

where F = i^K'^g^ is the decay rate of the vibrational amplitude. Here, the number of 
quanta in the bath mode n^, the coupling function K^, and the density of bath states g^^ are 
evaluated at a;^ = 2tc;, i.e., at the double frequency of the selected oscillator. 

Rewritten in the basis of eigenstates \n) of the unperturbed oscillator this master equation 
contains only linear combinations of such terms as (Tm,n — c Cm+2,n+2, and (Tm-2,n-2- 
It ensures, in effect, a possibility to distinguish even and odd initial states of the system. 
The odd excited state |1) cannot relax to the ground state |0), but the even excited state 
1 2) can. 

The evolution of the system after preparation under different initial conditions was sim- 
ulated numerically. The equations of motion of the density matrix elements are integrated 
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using a fourth order Runge-Kutta algorithm with stepsize control. To make the set of 
differential equations a finite one we restrict the number of levels by m, n < 20. 

Investigating the evolution of the initially coherent state of such a system we have con- 
centrated our attention on some detail, which was not seen before. One of the representative 
examples is the behavior of the initially coherent state of the system. For comparison with 
the usual behavior we show the time dependence of the mean value of the coordinate in Fig. 
4a. The mean value in the usual system decreases with a constant rate. The same initial 
value of the system coupled to a bath with cu^ — 2a; shows a fast decrement in the first stage 
and almost no decrement afterwards. 

The question about the influence of the type of the bath on the evolution of the superpo- 
sitional state is of special importance. At a temperature of fc^T = 2^a;/ln3, corresponding 
to n(2a;) = 0.5, we have simulated the evolution of the same superpositional states twice, 
for CO ^ — CO and cu^ — 2u!. The coordinate representation of the wave packets is presented in 
Fig. 4b and 4c. The same value of the relaxation 7 = F yields, however, different results. In 
the system coupled to the usual bath the quantum interference disappears already during 
the first period, while the amphtude decreases only slightly. In the other case, the system 
coupled to the bath with cu^ — 2a; leaves the interference almost unchanged, although a fast 
reduction of the amplitude occurs. 

Therefore, the second system partially conserves the quantum superpositional state. The 
experimental investigation of systems displaying such properties [2] is necessary, both for 
extracting typical parameter ranges for theoretical models and for practical applications like 
quantum computation and quantum cryptography. 

The described example leads us to the conclusion that the second type of environment 
is more preferable. 
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V. CONCLUSIONS 



Starting from von Neumann's equation for a vibrational oscillator interacting with the en- 
vironment modeled by a set of independent harmonic oscillators we derived a non-Markovian 
master equation, which has been solved analytically. For two types of the bath with maxi- 
mum of the spectral density near the system frequency and near twice the system frequency 
and for two different initial states, namely a coherent state and a superposition of coherent 
states, the wave packet dynamics in coordinate representation have been analyzed. It has 
been shown that wave packet dynamics demonstrate either "classical squeezing" and a de- 
crease of the effective vibrational oscillator frequency due to the phase-dependent interaction 
with the bath or a time-dependent relaxation rate, distinct for even and odd states, and par- 
tial conservation of quantum superposition due to the quadratic interaction with the bath. 
The decoherence also shows differences compared to the usual damping processes adopting 
RWA and to the description using the pointer basis for decoherence processes. We conclude 
that there is no permanent pointer basis for the decay. There are two universal stages of re- 
laxation: the coherence time scale stage and the Markovian stage of relaxation, both having 
different pointer bases. We believe that the proposed method can be applied for other initial 
states and different couplings with the environment in real existing quantum systems, which 
is important in the light of recent achievements in single molecule spectroscopy, trapped ion 
states engineering, and quantum computation. 
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FIGURES 

FIG. 1. Dynamics of the wave packet P{Q,t) (left plot) and its variance y(rigtit plot) for 
7 = O.lo;, = 3huj, Qo = 4. 

FIG. 2. Time evolution of the superposition of coherent states for a; = 1, 7 = O.Olo;, n = 0, 
a = 2, = f. 

FIG. 3. Time dependence of Pint {Q,t) for a; = 1, 7 = 0.25a;, n = 0.4, a = 2, = 0, Qo = 0. 
The solid line represents the case (31), the boxes represent RWA, the bullets represent earliest-time 
analysis (30) . 

FIG. 4. Influence of different baths, (a) Evolution of the coherent state: Mean value of the 
coordinate, for Qo = —2.2. Diamonds: system coupled to the bath with = u, j = 0.15a;. Solid 
line: system coupled to the bath with a;^ = 2a;, F = 0.5a;. (b) Time evolution of the superposition 
of coherent states with initial separation 2Qo = 8 for o;^ = a; = 1, 7 = 0.005a;, n = 1.36. (c) Same 
as (b), but for a;^ = 2a; with F = 0.005a;, = 0.5. 
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